%% Referentieimplementatie taylorsimple

function exp = taylorsimple(A)
exp=zeros(size(A));
add = eye(size(A));
% Tolerantie in taylorexpansie
tol=10^-10;
for i = 1:1000
    exp = exp+add;
    add = add*A/i;
    % Aanpassing om de relatieve fout te nemen
    if(norm(add/exp)<=tol)
        break
    end
end